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10 . In this article, we investigate the contribution of the high twist Feynman diagrams to the large-pr 

' single pseudoscalar and vector mesons inclusive production cross section in two-photon collisions 
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and we present the general formulae for the high and leading twist differential cross sections. The 
pion wave function where two non-trivial Gegenbauer coefficients 02 and 04 have been extracted 



QhI from the CLEO data, Braun-Filyanov pion wave function, the asymptotic and the Chernyak- 

Zhitnitsky wave functions are used in the calculations. For p-meson we used Ball-Braun wave 
function. The results of all the calculations reveal that the high twist cross sections, the ratio R, 
^ ! the dependence transverse momentum pT and the rapidity y of meson in the ^cleo{x, Q^) wave 

I function case is very close to the ^asyix) asymptotic wave function case. It is shown that the high 

twist contribution to the cross section depends on the choice of the meson wave functions. 
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I. INTRODUCTION 



During the last few years, a great deal of progress has been made in the investigation of 
the properties of hadronic wave functions[l-12]. The notion of distribution amplitudes refers 
to momentum fraction distributions of partons in meson in particular Fock state with fixed 
number of compenents. For the minimal number of constituents, the distribution amplitude 
$ is related to the Bethe-Salpeter wave function $^5 by 



The standard approach to distribution amplitudes, which is due to Brodsky and Lep- 
age[13], considers the hadron's parton decomposition in the infinite momentum frame. A 
conceptually different, but mathematically equivalent formalizm is the light-cone quantiza- 
tion[14]. Either way, power suppressed contributions to exclusive processes in QCD, which 
are commonly refered to as higher twist corections. The meson wave functions (also called 
distribution amplitudes -DA) [1] play a key role in the hard-scattering QCD processes be- 
cause they encapsulate the essential nonperturbative features of the meson's internal struc- 
ture in terms of the parton's longitudinal momentum fractions Xi. Meson wave functions 
have been extensively studied by using QCD sum rules. The original suggestion by Chernyak 
and Zhitnitsky of a "double-humped" wave function of the pion at a low scale, far from the 
asymptotic form, was based on an extraction of the first few moments from a standard QCD 
sum rule approach [5], which has been criticized and revised in Refs.[6,7]. Subsequently, 
a number of authors have proposed and studied the modified versions of meson [7,8] and 
baryon wave functions [9,10]. Additional arguments in favour of a form of the pion wave 
functions close to the asymptotic one have come from the analysis of the transition form 
factor 77* —>■ 7r° [12]. The measurements of this form factor by the CLEO collaboration are 
consistent with a near-asymptotic form of the wave function[15]. In [16], the leading- twist 
wave function of the pion at a low normalization point is calculated in the effective low- 
energy theory derived from the instanton vacuum. These results for the pion wave function 
at the low normalization point are close to the asymptotic form and consistent with the 
CLEO measurements. The authors have obtained a shape substantially different from the 
Chernyak-Zhitnitsky one because they have chosen a significantly smaller value of the second 
moment, and, more importantly, they have taken all the moments of the wave function into 
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account. Their results support the conclusions reached previously in Refs.[6,7]. The QCD 
factorization theorems predict that the hadron-hadron cross section can be obtained by the 
convolution of parton distribution functions and a cross section of the corresponding hard 
scattering subprocess. The parton distributions are nonperturbative, process-independent 
quantities, which are specific to any given hadron. The hard scattering cross sections are 
independent of all long distance effects and can be found by means of pQCD. In the frame- 
work of pQCD, the higher order corrections to the hard scattering, and therefore to the 
hadron-hadron process cross sections, have been calculated [17]. These corrections are large 
and change the leading order results considerably. Other corrections to the hadron-hadron 
process cross sections and its different characteristics come from the higher twist (HT) terms. 
By taking these points into account, it may be asserted that the analysis of the higher twist 
effects on the dependence of the meson wave function in single pseudoscalar and vector meson 
production at photon-photon collisions are significant in both theoretical and experimental 
studies. Much effort has recently been devoted to the study of exclusive processes involv- 
ing large transverse momenta within the context of perturbative quantum chromodynamics 
(QCD). Here, as in other applications of perturbative QCD, photon-induced reactions play 
an important role. In [18] have been studied photon-photon anihilation into two mesons at 
large center-of-mass angles Oc.m.- The higher-twist contributions to high-pr inclusive meson 
production in two-photon collisions, a single meson inclusive photoproduction and jet photo- 
production cross sections were studied by various authors [19-21]. As experiments examining 
high-pT particle production in two-photon collisions are improved, it becomes important to 
reassess the various contributions which arise in quantum chromodynamics. Predictins for 
the higher-twist contributions, orginally obtained in Ref.22, may now be refined using the 
exclusive-process QCD formalism devoloped in [23]. Another important aspect of this study 
is the choice of the meson model wave functions. In this respect, the contribution of the 
high twist Feynman diagrams to a single meson production cross section in photon-photon 
collisions has been computed by using various meson wave functions. Also, the leading 
and high twist contributions have been estimated and compared to each other. Within this 
context, this paper is organized as follows: in section [Tll we provide some formulae for the 
calculation of the contribution of the high twist diagrams. In section llllt we provide the 
formulae for the calculation of the contribution of the leading twist diagrams and in section 
IIVI we present the numerical results for the cross section and discuss the dependence of the 
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cross section on the meson wave functions. We state our conclusions in section |Vl 



II. CONTRIBUTION OF THE HIGH TWIST DIAGRAMS 

The high twist Feynman diagrams, which describe the subprocess 7g —>■ Mq contributes 
to 77 — s> MX for the meson production in the photon-photon colhsion are shown in Fig. 1(a). 
The amphtude for this subprocess can be found by means of the Brodsky-Lepage formula 
[24] 

M{s,i) = / dxi / dx26{l - Xi - X2)^M{xi,X2,Q'^)TH{s,i;xi,X2). (2.1) 
Jo Jo 

In Eq.(2.1), Th is the sum of the graphs contributing to the hard-scattering part of the 
subprocess. The hard-scattering part for the subprocess under consideration is 7g — > Mq, in 
which the observed meson is made directly. The hard-scattering amplitude T/y(s, t; xi, X2) 
depends on a process and can be obtained in the framework of pQCD, whereas the wave 
function $m(xi, X2, Q"^) describes all the non-perturbative and process- independent effects of 
hadronic binding. The hadron wave function gives the amplitude for finding partons (quarks, 
gluons) carrying the longitudinal fractional momenta x = (xi, X2, ■■■■Xn) and virtualness up to 

within the hadron and, in general, includes all Fock states with quantum numbers of the 
hadron. But only the lowest Fock state (gig2-for mesons, uud-ioi proton, etc.) contributes 
to the leading scahng behavior, other Fock state contributions are suppressed by powers 
of 1/Q^. In our work, we have restricted ourselves to considering the lowest Fock state 
for a meson. Then x = Xi,X2 and Xi + X2 = 1. This approach can be applied not only 
to the investigation of exclusive processes [25], but also to the calculation of higher twist 
corrections to some inclusive processes such as large-p^ dilepton production [26], two-jet -|- 
meson production in the electron-positron annihilation [27], ets. The qiq2 spin state used in 
computing Th may be written in the form 



E Usi[XiPm)Vs2[X2Pm) j^s 

V^lV^ ' '''' 



v/2 ' ^' 



^, PL heUntyO, (2.2) 
^ Pt heHcity ± 1, 

where e± = T(l/y2)(0, 1, ±i, 0) in a frame with (pm) 1,2 = and the N^^^^ project out 
a state of spins s, and pm is the four-momentum of the final meson. In our calculation, 
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we have neglected the meson mass. Turning to extracting the contributions of the high 
twist subprocesses, there are many kinds of leading twist subprocesses in 77 collisions as the 
background of the high twist subprocess 75 — > Mg, such as 7 + 7 q + q. The contributions 
from these leading twist subprocesses strongly depend on some phenomenological factors, 
for example, quark and gluon distribution functions in meson and fragmentation functions 
of various constituents etc. Most of these factors have not been well determined, neither 
theoretically nor experimentally. Thus they cause very large uncertainty in the computation 
of the cross section of process 77 — > MX. In general, the magnitude of this uncertainty is 
much larger than the sum of all the high twist contributions, so it is very difficult to extract 
the high twist contributions. 

The Mandelstam invariant variables for subprocesses jq — > Mq are defined as 



In our calculation, we have also neglected the quark masses. We have aimed to calculate 
the single meson production cross section and to fix the differences due to the use of various 
meson wave functions. We have used six different functions: the asymptotic wave function 
ASY, the Chernyak-Zhitnitsky [2,5] and the wave function in which two non-trivial Gegen- 
bauer coefficients 02 and 04 have been extracted from the CLEO data on the 77* — > 7r° 
transition form factor [28] and Braun-Filyanov wave function [7] . In ref.[28], the authors 
have used the QCD light-cone sum rules approach and have included into their analysis 
the NLO perturbative and twist-four corrections. They found that in the model with two 
nonasymptotic terms, at the scale /Iq — 2.AGeV. For p-meson we used Ball-Braun wave 
function [29]. In order to proceed to numerical calculations we have to use the explicit 
expressions for the p meson wave functions. They are defined by means of the following 
formulas: For the longitudinally polarized p meson, Pl = Ph=o 



(2.3) 



< I d{z)^^u{-z) I pl{p) >= Hp, f ciee^^(^^Vi(0, 



(2.4) 



J-1 



for the transversely polarized p meson, pT = Ph=±i 



< 



diz)a,Ju{-z) I pt{p) >= /J(e>. - e^p,) f dCe'^^'^'Um, (2-5) 



J-1 



fpi fp 2ire the dimensional constants which determine the values of the wave functions at 



the origin, and 




5 



In Eqs.(2.4), (2.5) and (2.6), d{z), u{z) and A^{a) are quark and gluon fields, is the 
polarization vector, and |pl(p)), \pt{p)) are the longitudinally and transversely polarized p 
meson states with the momentum p. 

^asyix) = VSf.xil - X), $^7^)(x) = V^/,^(^)x(l - x) 

^cz{x, fll) = 5^asy{2x - if, *^(^)(x, ^l) - ^t^T){x) [a + b{2x - If] , 

^BPix, l4) = ^asyix)[l + 0.66(5(2x - 1)' - 1) + 0.4687(21(2x - 1)^ - U{2x - if + 1)], 

<^CLEo{x,lil) = <^asy{x)[l + 0.285{5{2x-lf-l)-0.263{21{2x-lf-U{2x-lf + l)], (2.7) 

where /^=0.0923 GeV, fp—O.lAl GeV, /J=0.16 GeV is the pion and p mesons decay con- 
stants, a = 0.7, b = 1.5 for both longitudinally and transversely polarized p mcson[29]. 
Here, we have denoted by x = Xi, the longitudinal fractional momentum carried by the 
quark within the meson. Then, X2 = 1 — x and — 0:2 = 2a; — 1. The pion and p meson 
wave function is symmetric under replacement Xi — X2 X2 — Xi. The values of the pion 
wave function moments < > are defined as 

<r>= f^diCMO (2.8) 

Here, $7r(0 is the model function without and ^ = Xi — X2- The pion wave function 
moments have been calculated by means of the QCD sum rules method by Chernyak and 
Zhitnitsky at the normalization point po — O.BGeV. They are equal to 

< e >/.o= 1, < f >/.o= 0.44, < t >/.o= 0.27 (2.9) 
The Chernyak- Zhitnitsky pion model wave function has the following moments 

< e >Mo= 1, < e >Mo= 0.43, < >^„= 0.24 (2.10) 

It is interesting to note that the corresponding moments of the asymptotic wave function 
differ considerably from those in Eqs.(2.9), (2.10) 

< e >/.o= 1, < f >Mo= 0.20, < t >Mo= 0.086 (2.11) 

This means that the realistic pion wave function is much wider than the asymptotic one 
[5,30]. The model functions can be written as 

^asy{x) = V3f-;,X{1 - X), 

6 



cf{2x-i) + :^cr{2x-i) 



Co^\2x - 1) + 0. 18(0.2) -C2/"(2a; - 1) 

3 



.3/2. 



^bf{x, fJ,l) = ^asy{x) [Co^\2x - 1) + 0.44C|/'(2x - 1) + 0.25Cf '(2x - 1) 
^CLEo{x,lil) = <^asy{x) fco^'(2x - 1) + O.lQCa'/' (2x - 1) - 0.14Cf '(2x - 1) 



C',/\2x - 1) = 1, Cn^x - 1) = ^(5(2a; - if - 1), 



■y3/2/ 



Cf '(2a; - 1) = — (21(2x - 1)^ - 14(2a; - 1)' + 1). 
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(2.12) 



It may be seen that the pion wave function extracted from the experimental data depends 
on the methods used and their accuracy. Although one may claim that the meson wave 
function is a process-independent quantity, describing the internal structure of the meson 
itself, the exploration of different exclusive processes with the same meson leads to a variety 
of wave functions. This means that the methods employed have shortcomings or do not 
encompass all the mechanisms important for a given process. Such a situation is pronounced 
in the case of the pion. It is known that the meson wave function (distribution amplitude- 
DA) can be expanded over the eigenfunctions of the one- loop Brodsky- Lepage equation, i.e., 
in terms of the Gegenbauer polynomials {Cn^{2x — 1)}, 



1+ an{Q')C'J'{2x-l) 



n=2A.. 



(2.13) 



The evolution of the wave function (DA) on the factorization scale is governed by the 
functions a„(Q'), 



dniQ ) = an()Uo) 



(2.14) 



J2^50 74 ^364 
(3o 81' (3o 405' 
In Eq.(2.14), {jn} are anomalous dimensions defined by the expression. 



7n = Cf 



(n + l)(n + 2) 



n+l -I 



j=2 



(2.15) 



The constants ari(/io) = o-n are input parameters that form the shape of the wave functions 
and which can be extracted from experimental data or obtained from the nonperturbative 
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QCD computations at the normalization point //q. The QCD couphng constant ^^(Q^) at 
the two-loop approximation is given by the expression 



4:71 



(2.16) 



Here, A is the QCD scale parameter, /3o and /3i are the QCD beta function one- and two-loop 
coefficients, respectively, 

2 19 
/5o = 11 - 2^/, A = 51 - —Uf. 

For completennes we give the sum rules for the p- meson wave functions moments {i) — 
J d^e^u, 1^) [29]. 

(/p^)'(/^)(r)±e-"^?/^' = ^ £' ds £ due-^/^\u{2u - !)"• 

1 + f^(6 - ^ + 2/n4 + /nt. + Inn + In^l) ] + \j ^G-) + 

Stt 3 /x"' M J n -I- 1 12m^ tt 

''-(n-l)(^gg)^ (2.17) 



81M4' 



^"^^'^^ 47r2(n + l)(n + 3)^^ TT 

' {^G-)^^Mn-l){^-qq?. (2-18) 



12M2 ' TT ' 81M4 
Here the vacuum condensates are equal to [31] 

(^G^) = (0.012 ± 0.006)GeT/^ 

TT 

iV^mf = 0.56(-0.25G'e\/)^ 

The higher-twist subprocess 7g ^ Mq contributes to 77 — > MX through the diagram 
of Fig. 1(a). We now incorporate the higher-twist (HT) subprocess 7g Mq into the full 
inclusive cross section. In this subprocess 75 — >■ Mq, photon and the meson may be viewed 
as an effective current striking the incoming quark line. With this in mind, we write the 
complete cross section in formal analogy with deep-inelastic scattering, 

E-^{jj ^ MX) ^ -V^ / dx5{s + i + u)sGg/^{x,-t)^{^q^ Mq) + {t^u), (2.19) 

dp TT _ Jq dt 
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Here Gqj^ is the per color distribution function for a quark in a photon. The subprocess 
cross section for tt, pi, and pr production 



^(7? ^ Mq) 
at 



(2.20) 



where 



D{s, u) = ueias ( ^ ) (| ' + ^^2as 



where Ql — s/2, Q\ — — ■u/2, represents the momentum squared carried by the hard gluon 
in Fig. 1(a), 61(62) is the charge of ?i(?2) ^'^'^ Cp — \- 

The Im factors reflect the exclusive form factor of the meson and are discussed thoroughly 

in [20], as is the motivation the arguments of ag and Im- Note that the relation between Im 
and the meson form factor completely fixes the normalization of the higher-twist subprocess. 
The full cross section for tt and pL production is given by 



da 



MX) 



s -\-u 



s + t 



' 1 1 



P u'^ 



+ 



In (2.21), the subprocess invariants are 



(2.21) 



xs. 



t = t, 
u — xu, 

t = -2^^R - ^f) = -mT\/se~y, 



(2.22) 



with xr = (x|. + ,x|)^/^. Here xp = 2{pm)\\/^/s and xt = 2{pm)±/\/s = 2pt/\/s specify the 
longitudinal and transverse momentum of the meson. In terms of these the rapidity of M 
is given by 

y = \\.i^R + ^f)/ {xr - xf)] 
where uit - is the transverse mass of meson, which is given by 



As seen from (2.20) the subprocess cross section for longitudinal pl production is very 
similiar to that for tt production, but the transverse pr subprocess cross section has a quite 
different form. We have extracted the following high twist subprocesses contributing to the 

two covariant cross sections in Eq.(2.19) 

7?i ^ (?i?2)?i ,7?2 ^ (?i?2)?2 (2-23) 

As seen from Eq.(2.21), at fixed pt, the cross section falls very slowly with s. Also, at fixed 
s, the cross section decreases as l/pf^, multiplied by a slowly varying logarithmic function 
which vanishes at the phase-spase boundary. Thus, the pr spectrum is fairly independent 
of s expect near the kinematic limit. 



III. CONTRIBUTION OF THE LEADING TWIST DIAGRAMS 



Regarding the high twist corrections to the meson production cross section, a comparison 
of our results with leading twist contributions is crucial. The contribution from the leading- 
twist subprocess 77 — > qq is shown in Fig.l(b). The corresponding inclusive cross section 
for production of a meson M is given by 



da 



77^ MX 



E f'p{s + i + u)sDi^{z,-iAjj^qq) (3.1) 
- Jo z at 



where 



t ^ u 
s — s, t — - u — — 
z z 



Here s, t, and u refer to the overall 77 — > MX reaction. (z, —t) represents the quark 
fragmentation function into a meson containing a quark of the same flavor. For ti^ produc- 
tion we assume = In the leading twist subprocess, meson is indirectly emitted 
from the quark with fractional momentum z. The 5 function may be expressed in terms of 
the parton kinematic variables, and the z integration may then be done. The final form for 
the leading-twist contribution to the large-p^ meson production cross section in the process 
77 MX is 

da 3 x-^ dz 



■'M 



E 



d^P 



— JO 
?,9 



1 

-I 

Z 



■K ^ Jo z^ dt 

qq) 



d^6{s + ^{t + u))sDf{z, -t)^{ll 



da 
di 



l^^E-D, (z) 



t u 

- + ^ 

u t 



(3.2) 
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where 

t + u 
z — 

s 

We should note that D[z, —i) / z behaves as 1/2;^ as ;s — >■ 0. For the kinematic range 
considered in our numerical calculations, D[z, —i)/z increases even more rapidly. We obtain 
of the final cross section, Eq.(3.2), are as follows: At fixed pr, the cross section decreases 
with s asymptotically as 1/s. At fixed s, the D{z, —t) function causes the cross section to 
decrease rapidly as pt increases towards the phase-spase boundary [z ^ \). As s increases, 
the phase-spase boundary moves to higher p^, and the pt distribution broadens. 



IV. NUMERICAL RESULTS AND DISCUSSION 



In this section, the numerical results for higher twist effects on the dependence of the 
chosen meson wave functions in the process 77 — > MX are discussed. We have calculated 
the dependence on the meson wave functions for the high twist contribution to the large-pT 
single pseudoscalar tt"*" and vector p^, mesons production cross section in the photon- 
photon collision . The vr^, p^, p^ cross sections are, of course, identical. In the calculations, 
the asymptotic ^asy) Chernyak-Zhitnitsky ^czi Braun-Filyanov wave function[7], and also, 
the pion wave function, from which two non-trivial Gegenbauer coefficients 02 and 04 have 
been extracted from the CLEO data on the 7r°7 transition form factor have been used[20]. 
For p-meson we used Ball-Braun wave function[29]. In the ref.[28], authors have used the 
QCD light-cone sum rules approach and included into their analysis the NLO perturba- 
tive and twist-four corrections. For the high twist subprocess, we take 7g Mq and we 
have extracted the following two high twist subprocess 7^1 — > (^'1^2)91) IQ.2 ~^ {0.10.2)0.2 
contributing to 77 — > MX cross sections. Inclusive meson photoproduction represents a 
significant test case in which higher-twist terms dominate those of leading twist in certain 
kinematic domains. For the dominant leading twist subprocess for the meson production, we 
take the photon-photon annihilation 77 qq, in which the M meson is indirectly emitted 
from the quark. As an example for the quark distribution function inside the photon has 
been used [32]. The quark fragmentation function has been taken from [33]. The other 
problems dealth with are the choice of the QCD scale parameter A and the number of the 
active quark fiavors Uf. The high twist subprocesses probe the meson wave functions over 
a large range of squared momentum transfer, carried by the gluon. Therefore, we take 
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Qi = s/2, Q2 = —u/2 which we have obtained directly from the high twist subprocesses 
diagrams. The same has been used as an argument of as{Q'^) in the calculation of 
each diagram. The results of our numerical calculations are plotted in Figs. 2-7. Figs. 2- 
3 show the dependence of the differential cross sections of the high twist and ratio 

R = Y?IJ as a function of the meson transverse momentum for six different meson 
wave functions. As shown in Fig.2, the high twist differential cross section is monotoni- 
cally decreasing with an increase in the transverse momentum of the meson. As seen from 
Fig.2, in all wave functions of the mesons, the dependencies of the high twist cross sections 
on the "Pt transverse momentum of the meson demonstrate the same behavior. Also, as 
seen from Fig.2 the leading twist cross section is 2-4 order suppress the high twist cross 
section in magnitude, on the dependence wave functions of meson, respectively. On the 
other hand, the higher twist corrections are very sensitive to the choice of the meson wave 
function. We should note that the magnitude of the high twist cross section in the pion 
wave function ^cleo{^-, Q^) case is very close to the asymptotic wave function (^asy{x) case. 
In Fig. 3, the ratio R = is plotted at ?/ = as a function of the meson trans- 

verse momentum "Pt for the different meson wave functions. First of all, it is seen that the 
values of R for fixed y and \fs depend on the choice of the meson wave function. Also, 
the distinction between R{^asy{^)) with R{^cleo{^-,Q'^))-, Ri^czix^Q"^)), R{^bf{x,Q'^)), 
(x, Q^)) and Q^)) have been calculated. We have found that the distinction 

i?($asy(x)) and Ri^cLEoi.x^Q'^)) is small, whereas a distinction between R{^asy{,x)) with 
i?($cz(a;, Q^)), R{^BFix,Q'^)), Q^)) and i?($py(x, Q^)) is significant. For exam- 

ple, in the case of y/s = 89GeV, y = 0, the distinction between R{^asy{x)) with i?($j(a;, Q^)) 
{i = CLEO,CZ,BF,BB{L),BB(T)) is shown in Tabled Thus, the distinction between 
R{^asy{x)) and R{^cLEo{X: Q"^)) is maximum at px = l^GeV/c, but the distinction between 
R{^asy{x)) with R{^bf{x,Q'^)), R{^p^{x,Q^)) R{^p^{x,Q^)) is maximum 

at pt = 2GeV/c and decreases with an increase in pt- Such a behavior of R may be ex- 
plained by reducing all moments of the meson model wave functions to those of ^asy{x) 
for high Q^. In Fig. 4, the ratio R = Y?IJ jYJ^ is plotted at pt = 7GeV/c as a function 
of the rapidity y of the meson for the different meson wave functions. As we are now in 
the high energy region, the change of the rapidity to determine these relations is given by 
—ln{^/s/pT) < y < ln{^/pT)- At ^/s = 89GeV and pr = 7GeV/c, the meson rapidity lies 
in the region —2.543 < y < 2.543. First of all, it is seen that the values of R for fixed pt and 
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^/s depend on the choice of the meson wave function. As shown in Fig. 4 in all wave functions 
of the mesons, the dependencies of the ratio R = /T,'^ of the rapidity y of the meson has 
a minimum approximately at one point y = 1.75. After this point, the ratio increases with in- 
creasing y. Also, the distinction between R{^asy{x)) with R{^cleo{x, Q^)), R{^cz{x, Q'^)), 
R{^bf{x, Q^)), R{^pj^{x, Q'^)) and R{^p^{x, Q"^)) have been calculated. We have found that 
the distinction R{(^asy{.x)) and Ri^^cLEoiXiQ"^)) is small, whereas a distinction between 
R{^asy{x)) With R{^cz{x,Q^)), (o^, Q') ) , R{^p,{x,Q^)) and i?($,^(a;,Q2)) is gig. 

nificant. For example, in the case of y/s = 89GeV^, px = 7GeV/c, the distinction between 
R{^asy{x)) with R{^i{x, g2)) = CLEO, CZ, BF, BB{L), BB{T)) is presented in Table[TTl 
Thus, the distinction between R{(^asy{,x)) and R{^cleo{.x,Q'^)) is maximum at ?/ = —2.25, 
but the distinction between R{^asy{,x)) with i?($cz(a;, Q^)), R{^BFix,Q'^)), (x, Q^)) 

and R{^pj,{x,Q'^)) is maximum at y = 1.75 and decreases with an increase in y. Such a 
behavior of R may be explained by reducing all moments of the pion model wave functions 
to those of ^asy{x) for high Q^. We have also carried out comparative calculations in the 
center-of-mass energy ^/s = 209GeV. Figs. 5-6 show the dependence of the differential cross 
sections of the high twist S^-^, and ratio R = T,^^ /T,^ as a function of the meson trans- 
verse momentum for six different meson wave functions. As shown in Fig. 5, the high 
twist differential cross section is monotonically decreasing with an increase in the transverse 
momentum of the meson. As seen from Fig. 5, in all wave functions of the mesons, the 
dependencies of the high twist cross sections on the px transverse momentum of the meson 
demonstrate the same behavior. Also, as seen from Fig.5 the leading twist cross section is 4-5 
order suppress the high twist cross section in magnitude, on the dependence wave functions 
of meson, respectively. Also here, as in Fig. 2 the higher twist corrections are very sensitive to 
the choice of the meson wave function. We should note that the magnitude of the high twist 
cross section in the pion wave function ^cleo{XiQ^) case is very close to the asymptotic 
wave function $asy( In Fig.6, the ratio R = Sf//E^f is plotted at y = func- 

tion of the meson transverse momentum pt for the different meson wave functions. First of 
all, it is seen that the values of R for fixed y and ^/s depend on the choice of the meson wave 
function. Also, the distinction between R{^asy{x)) with R{^cleo{x,Q^)), R{^cz{x,Q'^)), 
R{^bf{x, Q"^)), -R($p^(x, Q^)) and -R($pj,(x, Q"^)) have been calculated. We have found that 
the distinction R{^asy{x)) and R{^cleo{x,Q^)) is small, whereas a distinction between 
R{^asy{x)) With R{^cz{x,Q^)), R{^ bf{x, Q^)) , R{^p,{x,Q^)) and R{^p^{x,Q^)) is signif- 
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icant. For example, in the case of = 2096*61^, y = 0, the distinction between R{^asy{x)) 
with R{^i{x,Q^)) {i = CLEO,CZ,BF,BB{L),BB{T)) is shown in Table [ml Thus, the 
distinction between R{^asyix)) with R{^cleo{.x, Q^)), R{^cz{,x, Q^)) and R{^bf{,x, Q"^)) is 
maximum at pt = 20GeV/c, but the distinction between R{^asy{LT){x)) with _R($p^(x, Q^)) 
and Q^)) is maximum at pt = lOOGeV/c and decreases with an increase in pt- In 

Fig.7, the ratio R = is plotted at pr = l7GeV/c function of the rapidity y of 

the meson for the different meson wave functions. At ^/s = 209GeV and pt = 17GeV/c, the 
meson rapidity lies in the region —2.509 < y < 2.509. First of all, it is seen that the values of 
R for fixed pt and ^/s depend on the choice of the meson wave function. Also, the distinction 
between R{^asy{x)) with R{^cLEoix,Q^)), i?($cz(x, Q^)), R{^bf{x,Q^)), R{^p,{x,Q^)) 
and -R($pj,(x, Q'^)) have been calculated. We have found that the distinction R{^asy{x)) and 
R{^CLEo{x,Q^)) is small, whereas a distinction between R{^asy{x)) with R{^cz{x,Q'^)), 
R{^bf{x,Q'^)), -R($p^ (x, Q^)) and R{^pj.{x,Q'^)) is significant. For example, in the case 
oi y/s = 209GeV^, pr = ViGeV/c, the distinction between R{^asy{x)) with Q^)) 
{i = CLEO,CZ,BF,BB{L),BB{T)) is shown in Table IVl Thus, the distinction between 
R{^asy{x)) and R{^cleo{x, Q"^)) is maximum at y = —2.25, but the distinction between 
R{^asy{x)) with R{^cz{x,Q^)), R{^bf{x,Q^)), R{^p^{x,Q^)) and R{^p^{x,Q^)) is maxi- 
mum at ?/ = 1.75 and decreases with an increase in y. Also, as seen from Fig.7 in all wave 
functions of the mesons, the dependencies of the ratio R = of the rapidity y of 

the meson has a minimum approximately at one point y = 1.75. After this point, the ratio 
increases with increasing y. As seen from calculations with increasing center-of-mass energy 
from y/s = S9GeV to y/s = 209GeV^, the distinction between R decreases for all meson wave 
functions. 

V. CONCLUDING REMARKS 

In this work, we have calculated the higher twist contribution to the large-pr meson 
production cross section to show the dependence on the chosen meson wave functions in 
the process 77 —>■ MX. In our calculations, we have used the asymptotic ^asy, Chernyak- 
Zhitnitsky ^cz, Braun-Filyanov ^bf wave functions and also, the pion wave function, in 
which the coefficients 02 and 04 have been extracted from the CLEO data on the tc^'j 
transition form factor used. For p-meson we used Ball-Braun wave function. For the high 
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twist subprocess, we have taken — > Mq. We have extracted the following two high twist 
subprocesses 7^1 — > {qiq-ijqi-i 7?2 ~^ {0.10.2)0.2 contributing to 77 MX cross sections. 
As the dominant leading twist subprocess for the meson production, we have taken the 
photon-photon annihilation 77 qq, where the M meson is indirectly emitted from the 
quark. The results of our numerical calculations have been plotted in Figs. 2-7. As shown in 
Figs. 2, 5 the high twist differential cross section monotonically decrease when the transverse 
momentum of the meson increases. As seen from Figs. 2, 5 in all wave functions of mesons, 
the dependencies of the high twist cross sections on the px transverse momentum of the 
meson demonstrate the same behavior. But, the higher twist corrections are very sensitive 
to the choice of the meson wave function. It should be noted that the magnitude of the high 
twist cross section for the pion wave function ^cleo{x,Q^) is very close to the asymptotic 
wave function ^asy{-i^)- 

In Figs. 3 and 6, the ratio R — jYJ^ has been plotted at y = as a func- 
tion of the meson transverse momentum, p^, for the different meson wave functions. 
It may be observed that the values of R for fixed y and ^/s depend on the choice of 
meson wave function. Within this context, we have also calculated the distinction be- 
tween R{^asy{x)) and i?($cLi?o(a;,g')), i?($cz(a:,Q2)), Q^))^ Q^))^ 
i?(<I>py(a;, Q^)) . We have ultimately found that the difference between i?(<l>asj/(a^)) and 
R{^CLEo{^-,Q^y) is small, whereas a distinction between R{^asy{^y) with R{^cz{^-,Q'^y)-, 
R{^bf{x,Q^)), R{^pj^{x,Q^)) and i?($pj,(x, Q^)) is significant. In Figs.4 and 7, the ratio 
it! = E^^/E^ has been plotted at pr — 7, llGeV/c as a function of the rapidity of the 
meson for the different meson wave functions. It may be observed that the values of R for 
fixed Pt and ^/s depend on the choice of meson wave function. Within this context, we have 
also calculated the distinction between R{^asy{x)) and R{(^cleo{x-,Q'^))-, Ri^cz{x,Q'^)), 
R{^bf{x,Q^)), R{^p^{x,Q^)), R{^p^{x,Q^)) . We have ultimately found that the dif- 
ference between R{^asy{x)) and R{^cleo{x,Q'^)) is small, whereas a distinction between 
R{^asy{x)) with R{^cz{x,Q^)), R{^bf{x, Q^)) , R{^p^{x,Q^)) and i?($p^(x, Q^)) is signif- 
icant. Our investigation enables us to conclude that the high twist meson production cross 
section in the photon-photon collisions depends on the form of the meson model wave func- 
tions and may be used for their study. Further investigations are needed in order to clarify 
the role of high twist effects in QCD. 
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VT,GeV/c 






R('S'bf(x,Q-')) 


R{'S>BBiL)(x,Q'^)) 


i?(*BB(T)(x,Q2)) 






K[<Pasy(x)) 


H(<Pasy(X)) 


K[<Pasy[x)) 


2 


0.232 


9.054 


23.999 


3.490 


4.048 


6 


1.492 


4.216 


5.872 


3.950 


2.861 


14 


1.773 


1.747 


1.216 


1.394 


1.537 


22 


0.9499 


0.7607 


0.68177 


0.85668 


0.80928 


28 


0.5568 


0.5023 


0.7200 


0.68469 


0.5859 


36 


0.9874 


0.79656 


0.6965 


0.8766 


0.8341 



TABLE I: The distinction between R{^asy(.x)) with R{^i{x,Q'^)) (i=CLEO, CZ, BF, BB(L), 
BB(T) ) at cm. energy = 89GeV. 



V 


R{^CLE0{X,Q'^)) 


R{^cz{x,Q'^)) 


R{^bf(x,Q'')) 


R{'i>BB{L){x,Q^)) 


R{<i>BB{T){x,Q^)) 






R{^asv{x)) 


R{^asv[x)) 


R[^asy[x)) 


-2.25 


1.545 


1.22 


0.745 


1.113 


1.146 


-1.5 


0.769 


0.594 


0.615 


0.773 


0.714 


0.25 


1.428 


4.183 


5.97 


2.399 


2.893 


1.0 


0.916 


5.169 


9.693 


2.773 


3.276 


1.75 


0.559 


6.936 


15.798 


3.186 


5.549 


2.25 


0.892 


4.922 


9.281 


2.744 


3.529 



TABLE II: The distinction between R{^asy{x)) with R{^i{x,Q'^)) (i=CLEO, CZ, BF, BB(L), 
BB(T)) at cm. energy ^ = 89GeF. 
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VT,GeV/c 






R('S'bf(x,Q-')) 


R{'S>BBiL)(x,Q'^)) 


i?(*BB(T)(x,Q2)) 


H['Pa3y(x)) 


H['Pa3y(X)) 


K(<Pasy(x)) 


K(<Pasy(X)) 


K[<Pasy(x)) 


20 


1.686 


2.760 


2.865 


1.914 


2.298 


35 


1.582 


1.507 


1.088 


1.294 


1.413 


50 


1.014 


0.835 


0.7266 


0.8962 


0.8577 


65 


0.614 


0.5534 


0.7367 


0.70467 


0.6017 


85 


1.011 


0.839 


0.736 


0.0.8977 


0.8586 


100 


1.574 


3.368 


4.118 


2.170 


2.637 



TABLE III: The distinction between R{^asy{x)) with R{^i{x,Q'^)) (i=CLEO, CZ, BF, BB(L), 
BB(T) ) at cm. energy = 209Gey. 



y 


R{^CLE0ix,Q'')) 
R{^asy{x)) 


R{^cz{x,Q^)} 

R{^asy{x)) 


R{^bf(x,Q'')) 

R{^asy{x)) 


R{^BB(L){X,Q'')) 
R{^asy{x)) 


R{^BB(T){X,Q'')) 
R{^asy{x)) 


-2.25 


1.571 


1.383 


0.933 


1.212 


1.289 


-1.5 


0.7519 


0.616 


0.673 


0.7689 


0.6978 


0.25 


1.435 


3.484 


4.626 


2.236 


2.754 


1.0 


1.039 


4.255 


7.252 


2.575 


3.149 


1.75 


0.7246 


5.887 


12.269 


2.999 


6.338 


2.25 


1.056 


4.003 


6.684 


2.520 


3.295 



TABLE IV: The distinction between R{^asy{x)) with R{^i{x,Q'^)) (i=CLEO, CZ, BF, BB(L), 
BB(T)) at cm. energy = 209Gel^. 
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FIG. 2: High twist meson M production cross sections as a function of the pr transverse momentum 
of the meson at the c.m.energy -s/s = 89 GeV. 
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FIG. 3: Ratio R = H^/Ttj^, where the leading and the high twist contributions are calculated for 
the meson rapidity y = at the cm. energy ^/s = 89 GeV, as a function of the meson transverse 
momentum, pT- 
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FIG. 4: Ratio R = /T,^ , where the leading and the high twist contributions are calculated 

for the meson transverse momentum px = IGeV/c at the cm. energy -y/s = 89 GeV , as a function 

of the y rapidity of the meson. 
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FIG. 5: High twist meson M production cross sections as a function of the px transverse momentum 



of the meson at the c.m.energy = 209 GeV . 
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FIG. 6: Ratio R = fEj^J, where the leading and the high twist contributions are calculated for 
the meson rapidity y = at the cm. energy ^/s = 209 GeV, as a function of the meson transverse 
momentum, px- 
10" 




FIG. 7: Ratio R = S^-^/S^, where the leading and the high twist contributions are calculated for 
the meson transverse momentum px = 17GeV/c at the cm. energy ^/s = 209 GeV, as a function 
of the y rapidity of the meson. 
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